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ABSTRACT 


Near infrared observations in the reflection nebulae NGC 7023, 
2023, and 2068 are used to study cloisters of young stars found 
associated with these nebulae. At least 30*60 % of these stars are 
p^e-main sequence objects, as indicated by their infrared 
excesses, hydrogen line emission, or irregular variability. The 
spatial distributions and observed luminosity functions of these 
young open clusters are derived, and the inferred mass function 
and star formation efficiencies are discussed. 





1. Observations made partially at the Mt. Wilson Observatory, Carnegie 
Institution of Washington, as part of a collaborative agreement between 
Carnegie Institution of Washington and California Institute of Technology. 
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L INTRODUCTION 

Near Infrared observations of dark clouds provide a unique approach to 
studying low mass star formation. Optical observations are hampered by the 
large extinctions associated with dark clouds, while radio continuum 
observations a'e sensitive only to more massive ionizing stars. Early near 
infrared searches for embedded stars found an obscured young cluster in the 
Ophiuchus dark cloud (Grasdalea Strom, and Strom 1973; Vrba et of. 1975), and 
searches were soon made in many other dark clouds (see Hyland 1961 for a 
review). Elias (1978a, 6,c) showed that many sources found in 2.2 /xm scans of 
dark clouds were background stars not associated with the regions, and that 
many stars associated with the clouds were pre-main sequence objects. The 
identification o? a substantial background component to 2.2 jum star counts in 
dark clouds led to an improved understanding of the near infrared extinction in 
these clouds. The presence of background stairs, and of stars not yet on the 
main sequence, also showed that extreme care must be taken in interpreting 
near infrared star formation studies in these regions. These effects have been 
considered in recent estimates of the star formation efficiency in the 
Chamaeleon dark cloud (Hyland, Jones, and Mitchell 1982). 

The observations discussed in this paper are of stars found in the regions 
near three reflection nebulae, NGC 7023, 2023, and 2066. These stars are shown 
to be members of young clusters associated with the reflection nebulae. The 
spatial distributions and luminosity functions of these stars are derived, and 
estimates of the mass functions and star formation efficiencies are made. The 
high success rate in finding young clusters of low mass stars is attributable to 
selecting reflection nebulae for these observations, since this guarantees the 
absence of high mass (0) stars which would ionize substantial H II regions, while 
identifying a region of star formation by the presence of a lower mass (B) early 
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type star still young enough to be physically associated with the cloud from 
which it formed. 

NGC 7023 is at the northern end of a 30' x 60' molecular cloud at a distance 
of 440 pc (Viotti 1969). The visual reflection nebula is illuminated by the Herbig 
Be star HD 200775, and the nebula is also associated with a small cluster of T 
Tauri stars (Weston 1953). NGC 2023 is found in the dark cloud L 1630 (Lynds 
1962), and is 20’ south oJ the H II region NGC 2024. The visual reflection nebula 
is illuminated by the B1.5V star HD 37903. NGC 2066, or M78, is also in L1630. 
NGC 2068 is 14* south of the active star formation region NGC 2071, and is also 
18' north of the Herbig-Haro object HH-24. The visual reflection nebula is 
illuminated by the northernmost of a pair of stars having V «- 10 mag (Elvius and 
Hall 1967). These stars, following the notation of Strom at al. (1975), will be 
referred to as HD 38563-N and HD 38563-S. NGC 2023 and 2053 are at a distance 
of 500 pc (Lee 1968). 

The observations reported here include scans at 2.2 yu.m, and occasionally 
1.65 pm, to search for stellar sources; near infrared photometry at 1.25, 1.65, 
and 2.2 pm for the brighter sources found at 2.2 pm; visual photometry for stars 
in two of the regions at 0.655 pm; visual and near infrared spectrophotometry 
for a few stars; and an image tube plate at 0.8 pm of one region. This diverse 
body of data allows an investigation of the colors, reddening, spectral features, 
and variability of the stars in these dark clouds. These observations lead to a 
better understanding of the nature of young stellar objects and low mass star 
formation oc curing in dark clouds, as well as the nature of the dark cloud 


extinction. 
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D. OBSERVATIONS 

Scant * A the reflection nebulae NGC 7023, NGC 2023, and NGC 2068 were 
made at 2.2 /im. and in some cases, 1.65 /un, using the 100 inch, 60 inch, and 24 
inch telescopes at Mt. Wilson Observatory. The areas scanned, wavelength, 
diaphragm, spacing between source and sky positions, and faintest magnitude 
detectable are listed in Table 1 for the three nebulae. Adjacent scans were 
spaced by half the diaphragm diameter. The scans were usually made in 
nonphotometric conditions, so that the inferred magnitudes are typically 
reliable to only ~ 0.3 mag, with worse uncertainties in some cases. 

A plate of the NGC 7023 region was obtained on the Yerkes Observatory 41 
inch reflector, at an effective wavelength of 0.6 /m i, using a KG-8 Alter and a S-25 
image intensifier tube. The 20 brightest stars from this plate were selected and 
searched for at infrared wavelengths. 

Infrared photometry at J (1.25 /im). H (1.65 /mi). and AT (2.2/xm), using the 
filters defined by Neugebauer at at. (1979), was obtained of many of the stars 
found in scans and on the 0.8 /mi plate. The photometry was calibrated relative 
to the standard stars given by Elias at at. (1982). 

Photometry at a wavelength of 0.655 /mi, using an r filter as defined by 
Thuan and Gunn (1976), was obtained for stars in NGC 2023 and NGC 2068 on 
1981 September 25. This photometry was derived from 11 direct SIT frames 
taken on the Palomar 60 inch telescope, using standard techniques described by 
Kent (1979). Nonphotometric sky conditions resulted in photometric 
uncertainties of ± 0.2 mag. 

A spectrum with a resolution of AX/ A = 1.7 x 10 -3 . from 2.15 to 2.22 /tm, was 
obtained of HD 200775 on 1980 October 9 with the Mt. Wilson Observatory 60 inch 
telescope. This spectrum was measured with a Fabry-Perot spectrometer, using 
techniques described in more detail by Sellgren at al. ( 1983a). 
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Spectra from 0.36 p m to 0.76 pm were obtained of 3 stars In NGC 3023 and 
NGC 3066 on 1961 August 30-31, using the Palomar 300 inch telescope and the 
double spectrograph blue camera described by Oke and Gunn (1963). 


0 • 


ID. RESULTS 
a. Aotonutry 

lb* observation* of atara in NGC 7023, 2023, and 2060 a r* praaantad in 
Table 2. The designation given in the Drat column will be the notation used to 
refer to the atara throughout this paper, although other designations are also 
given in Table 2 when available. 

Table 2 contains 30, 14, and 56 stars found in NGC 7023, 2023, and 2068 
respectively at infrared wavelengths, as well as some stars found only at visual 
wavelengths. Of these stars, respectively 11, 8, and 15 are brighter than K * 
10.5 mag in scan areas typically 6-8’ on a side. 

The background star count model of Jones • t aL (1981) can be used to 
estimate the expected number of background stars in these scan areas. This 
model agrees very well with observations of background stars in many regions. 
In particular scans of a 2.4 square degree region near NGC 7023 (Elias 1983, 
private communication) find 19.2 ± 2.8 stars per square degree brighter than K 
= 7.5 mag, while the model predicts 15.4 stars per square degree brighter than 
K = 7.5 mag, in good agr sement with the observations. The background star 
count model predicts that at the galactic latitudes of NGC 7023, 2023, and 2068, 
14.2*, 16.6*, and 14.5*, one expects 120, 75, and 86 background stars per square 
degree with K * 10 mag, respectively. Similarly one expects 230, 150, and 170 
background stars per square degree, respectively, with K « 11 mag. Thus in the 
areas scanned only 2 or less background stars with K * 10.5 mag are expected. 
Many more stars than this were found in each nebula, so that 75-90% of the 
stars found in near infrared scans are stars pnysically associated with the 
reflection nebulae. This confirms the suggestions of other observers of NGC 7023 
(Weston 1953; Strom si at. 1972), NGC 2023 (Strom si at. 1975), and NGC 2068 
(Herbig and Kuhi 1963; Strom si at. 1975; Strom. Strom and Vrba 1976) that 
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thesa regions contain clusters of young atara. Furthar arguments for a largo 
traction of atara phyiically aaaociatad with tha clouds art tha olaar spatial 
dusts ring of tha stars and tha idantlficatlon of many stars as pre-main saquanca 
stars, as dlscussad in mors detail in saotlons IVb and IVe . 

Photometry at /, H, and K was obtalnad for aoma of thosa stars, and la 
glvan in Table 2. All /, H, and K photometry for each star was obtalnad on tha 
same night, so that possible variability of a star should not affect Its infrared 
colors. The J-H and H-K stellar colors are plotted in Fig. 1 for each nebula. 
In these J -H vs. H-K diagrams, the reddening line (cf. section IVa) is drawn 
through a point corresponding to the average color of a foreground or 
background star, which model predictions at 2.2 /im indicate is most likely a K 
giant. Tha reddening line for an early type main saquanca star is also shown. 

In NGC 7023 (Fig. la) most of the stars appear to have colors appropriate to 
a normal reddened star. Three stars, however, fall substantially below the 
reddening line for early type stars, implying that these stars have an infrared 
excess. The remaining stars are, by their colors, consistent with either being 
background giants or stars without strong infrared excesses associated with the 
cloud. In NGC 2023 (Fig. lb), almost half the stars observed tall below the 
reddening line for early type stars, and therefore have infrared excesses. The 
color-color diagram for NGC 2066 is shown in Fig. lc. Of the stars in NGC 2066 
whose photometry is reported in Table 2, stars 5, 16, 19, and 30 have infrared 
excesses. For NGC 2066 the photometry of this paper was supplemented by the 
observations of Strom. Strom, and Vrba (1976), transformed to the photometric 
system of Neugebauer *t at. (1979) using the relation of Frogel cf at. (1978). Of 
the stars listed in Table 2 for which the colors at Strom ef at. have been 
adopted, stars 4 and 2 have infrared excesses, while stars 20 and 53 have no 
evidence for infrared excesses. All stars in the three nebulae with identified 



infrared axcauu have this noted in Table 2. 


The r photometry in NGC 2023 and 2068 is also presented in Table 2. 
Typically t-K ~ 6 mag In the two nebulae. The stars detected at only r or K, 
but net both, have upper or lower limits onr-JT generally consistent with r -K ~ 

6 mag also. 

6. Spac.traphotomatry 

A spectrum of HD 200775 was obtained from 2.15 to 2.22 /sm with AX/X = 1.7 
x 10"®. B y was seen in emission from the star, with a strength of 2.3 x 10" 1 * erg 
cm" 1 s' 1 , and was unresolved with an instrumental resolution of 500 1cm s" 1 . A 3 
a upper limit of 9 x 10~ 18 erg cm"* s" 1 was placed on other possible features in 
the spectrum. 

Visual spectrophotometry was obtained for two stars in NGC 2023, stars H 
and C, and for we star in NGC 2068. star 20. The continue of these stars are 
extremely red. with B-V = 1.6, 1.8, and 2.2 mag respectively for 2023-H, 2023-C, 
and 2068-20. In 2023-H and 2023-C, H a is seen in emission from the star, while 
in 2068-20 both H a and H ft are seen in emission. H a emission was previously 
identified in 2068-20 by Herbig and Kuhi (1963), who refer to this star as IicHa 
304. 

The hydrogen line emission in HD 200775, 2023-H, 2023-C. and 2068-20 must 
come from dense circumstellar shells or stellar winds. This is because the 
observed B y and H a fluxes would predict optically thin free-free emission at 
radio wavelengths, far in excess of the 6 cm upper limits observed by Sellgren 
at ol. (19836) for these stars. 

c. \hrlability 

The observations reported here also include measurements at K. and 
sometimes H, made at different times, which can be used to investigate the 
variability of the stars observed in the three reflection nebulae. Both scan 
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magnitudes tod photometric magnitudes obtained on different data* ara 
presented in Tabla 3. at wall aa maaauramanta by othar obaarvars whan 
ayailabla. Only atari auapactad of variability ara includad in this tabla. 

In comparison! batwaan different seta of photometry, with uncertainties ~ 
0.05 mag. difference a 0.15 mag ware considered evidence of variability. Whan 
comparisons includad Kan magnitudes, which ara uncertain by fc 0.3 mag. 
differences n 0.75 mag ware taken aa evidence of variability. Only differences a 
0.75 mag ware considered in establishing variability based on comparisons 
between the measurements presented in this paper and those of other 
observers, due to the possibility of unknown systematic differences between the 
various sets of data. 

In NGC 7023, 6 stars are definitely identified as variable stars, based on 
photometry, while another 3 stars may be variable, based on scan magnitudes. 
Only the stars whose photometric magnitudes varied are identified aa variable in 
Table 2a . These 6 stars varied by amounts ranging from 0. 17 mag to 0.59 mag at 
K. The shortest time scale variability was seen in star D, where two photometric 
measurements on 1961 July 26. separated by 5 hours, differed by 0.15 ± 0. J4 
mag at K. HD 200775 (star A) was extensively observed at /, H, and K. 
Although this star is a Herbig Be star, stars which as a class are frequently 
variable, HD 200775 Itself shows no evidence for infrared variability on time 
scales of hours, days, or years. 

There is Uttle evidence of variability * 0.75 mag for the stars In NGC 2023. 
Star F appeared to vary by 1.4 mag at K from Intercomparisons of different 
scans and photometry, but this is probably due to confusion (see footnote to 
Table 2b). 

In NGC 2066 13 stars may have varied by > 0.75 mag. Evidence for 
variability in stars 13 and 15 depends on single bright K measurements on the 
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Ht Wilson 24 inch telescope, near tha scan dataction limit. Tha variation In 
atari 26a and 26b la likaly due to confusion. Tha remaining 0 atari have baan 
identified aa variable in Table 2e. Star 14 varied 0.9 mag at H. baaed on 
me aiura manta of other observer! (Strom at at. 1975; Strom, Strom, and Vrba 
1976). The other 8 stars varied by amounts ranging from 1.0 mag to 1.5 mag at 
K. 

It is important to realize that while the fractions of stars found to vary in 
NGC 2023 and 2068 are less than in NGC 7023, the variability identified in NGC 
7023 would have been missed if two sets of photometry separated by two years 
had not been available. Only one set of photometry exists in the other two 
nebulae, so that only variability fe 0.75 mag can be detected by comparison with 
scan magnitudes. The largest variation in the NGC 7023 photometry is smaller 
than this. This suggests that additional photometry in NGC 2023 and 2068 may 
uncover many more variable stars. 
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lY.HSCUBBOM 

a. infrarad axtinction in refaction nabuloa 

Tte reddening line In Fig. 1, E{J-H)/E{H-K) ■ 1.90, la derived from tha 
value of Jonaa and Hyland (1960), maaaurad on tha AA0 photometric system, by 
transforming It (Ellas at at. 1963) to the photometric Altars used hero. A 
slightly different value, E{J-H)/E(H-K) » 1.6, has bean obtained by Ellas 
(1976c). Use of this reddening law instead would hare no affect on tha numbers 
of pro-main sequence stars found or on other conclusions of this paper. 

Jonas and Hyland (i960) have argued that the extinction curve in the near 
Infrared, at least at J, H , K, and perhaps L (3.4 /an). Is identical, within the 
uncertainties. In the general Interstellar medium and in dense molecular clouds 
where large values of R * Ay/E(B-Y) have been reported. As discussed by 
Hyland (1961), the upper envelope of a J-H vs. H-K diagram for the stars 
observed towards a dark cloud Is composed of background late type giants and 
stars without infrared excesses associated with the cloud, and therefore defines 
the near infrared extinction curve In the dark cloud. As can be seen from the 
good agreement between the reddening line of Jones and Hyland (1960) 
indicated in Fig. 1 and tha observed upper envelopes to the color-color 
diagrams, the observed near infrared extinction curves in these three reflection 
nebulae support the hypothesis that a single near infrared extinction curve 
applies in a wide variety of physical environments. 

While the near infrared extinction curve, as defined by E(J -H)/E[H-K), is 
seen to be constant, the absolute value of the extinction Is not known. Jonas and 
Hyland (1960) fit a power law to their infrared extinction curve at /. H, K. and 
L, and find that at these wavelengths that tha extinction is consistent with a 
XT** power law. This neglects the possibility of a grey component to the near 
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infrared extinction. Without more Information mi the absolute value at the mat 
infrared extinction, it seems reasonable to adopt on extinction proportional to 
X**, which Implies that S(J-H)/Ag * 2.06 and S(H~K) ! Ag • 1.06. Thsss 
values srs adopted throughout this paper. 

i. The fraction af start which art prt-matn nqutnct 

To determine the traction at the total number of stars which are pre-main 
sequence (PICS), one must first restrict the sample so that It Is complete and 
reasonably uniform. For this reason the only stars considered are those found 
In 2.8 pm scans in 1879 of NGC 7083, with a detection limit ot 11.0 mag; 8.8 pm 
scans in 1061 January of NGC 8083, with a detection limit ot ll.l mag; and 2.2 
pm scans in 1861 September and 1961 November of NGC 8066, with detection 
limits of 11.9 • 18.6 mag. These samples are estimated to be complete to K * 
10.S mag for NGC 7083 and 2023, and to 11.89 mag for NGC 2068. Start below 
the scan completeness limit in each case are excluded from the sample. This 
provides uniform selections of stare in each nebula found from a set of scans 
with a single completeness limit, over a wide area. The small fraction ot 
background stars, 10 - 89X, will be neglected in determining the characteristics 
of the stellar dusters associated with the reflection nebulae. The scan areas for 
these samples are 0.98, 1.39, and 1.86 pc* respectively for NGC 7023, 2023, and 
8068. The stars included in the sample are all stars in Table 8 brighter then the 
scan completeness limit, except for stars 40 and 48 in NGC 7023, stars L and N in 
NGC 2023, and stars 1. 2, 3, 4. 9, 6, and 96 in NGC 2068. Note that star D in NGC 
7023 Is Included in the sample, since it was brighter than the scan completeness 
limit when first discovered (Table 3a). 

Some of the stars which are listed in Table 2a are else identified with 
variable or suspected variable stars found by Rosino and Romano (1962). When 
suitable observations exist at both wavelengths, all of the stars known to vary in 
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NOC 7033 are seen to be variable at both visual and naa r infrared wavelengths. 
The study of Roatno and Romano (1003), which is more extensive la tt 'i Uma 
coverage. shows that tb» stars in NOC 7033 wary irregularly with Uma, which Is 
characterisUc of nebular variables such as T Taurl stars, fhe Uma coverage of 
tha observations roportad in Tabla 3 is not ooiLplata enough to aatabliah tha 
nature of tha variability of tha start observed In tha thraa reflection nabulaa. 
Tha good corraspondanco batwaan variability at K and IdsntificaUon as an 
irragular variable at visual wavalangths In NGC 7033, howsvar, suggests that it Is 
a reasonable assumption that infrared variability Indicates that a star is a PMS 
star physically associated with tha cloud. 

The variety of observations described above provide an approach to 
estimating the number of young PMS stars In the raflaction nabula fields. One 
method of identifying a PMS star is by its variability. The presence of an 
infrared excess or Ha emission, both of which are indicators of eircumstellar 
shells which surround PMS stars (Strom, Strom, and Grasdalen 1976), can also 
be used to identify PMS stars. These three characterisUcs are noted in Table 3 
for each star that is known to show them. 

In the sample wb'xh has been restricted so as to be complete, 6 out of 11 
stars in NGC 7023 are Identified as PMS stars, using the criteria of variability, 
infrared excesses, and h/drogen Una emission. Similarly 5 out of 3 stars in NGC 
2033, and 9 out of 28 stars Li NGC 3068, are PMS stars. Many of the irregular 
variables in NGC 7033 show no evidence of an infrared excess in the color-color 
diagrams, but the observations of this paper are not sufficient to determine the 
variability of moet stars in NGC 3033 and 3088. Thus these numbers set lower 
limits to the fraction ol the detected stars which are PMS stars. 

All the stars an NGC 3033 and 3088 cannot be PMS. since star A in NGC 3033 
is classified B1.5V (Abt and * evato 1977), and stars 31 and 36 in NGC 3068 are 
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clas sited ~ B3V (Sellgren at al. 19636). Thasa stars, in addition to star 33 in 
NGC 2068 (B2H-U1; Strom at at. 1975) appear to ba normal stars with no 
evidence of variability, H a emission, or infrared excesses. The earliest type 
star in NGC 7023, HD 200775, is a PMS star (Strom at at. 1972), with an infrared 
excess and H a emission, although it is not photometrically variable (section 
Hie; Breger 1974). The remaining stars in NGC 7023 are either PUS or 
unclassified, so it is possible that no stars in NGC 7023 have reached the main 
sequence yet. 

c. Jhs luminosity function and spatial distribution of ths stars 

The spatial distibution of stars in the three nebulae is shown in Fig 2. The 
surface density of stars can be seen to be highest near the cluster centers, 
which in NGC 2068 and NGC 7023 is consistent with the positions of the visual 
ill uminatin g stars of the reflection nebulae, star 33 and star A respectively. In 
NGC 2023 the visual illuminating star, star A does not appear to be at the center 
of the cluster, so the NGC 2023 cluster center is assumed to be the mean 
position of the stars, 60” east and 67” south of star A The radial distribution of 
stars firom the cluster center for each nebula appears to follow a power law of 
the form N(R) ~ R*. where N(R) is the total number of stars inside of a circle 
of projected radius R. Power laws were found tor NGC 7023, 2023, and 2066 of x 
= 0.94 ± 0.26, z = 0.63 t 0.33, and z = 0.B9 ± 0.10, respectively, using only the 
stars in the samples which have been limited so as to be complete. These are all 
consistent with z = 1, or N(R) ~ R. This distribution corresponds to a surface 
density of stars a(R) ~ R~ x , if a is uniform. If the clusters are spherically 
symmetric, the observed distribution N{R) also corresponds to a volume 
density of stars n(r) ~r~* where r is the true radius from the cluster center. 

The distribution in apparent magnitude can also be determined from the 
observations of the three reflection nebulae. The cumulative number of stars 
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N{K) m. observed K magnitude la shown in Tig. 3 for NGC 3068. where the 
largest number of start were found. This figure Includes data from all the K 
nans, up to a K magnitude 0.5 mag brighter than each scan detection limit. 
Mated to account for the different surface densities of stars sampled by 
different scans, and weighted by the uncertainties due to counting statistics In 
each magnitude bin. The assumption is made that N(K) Is Independent of 
position In the nebula. The observed distribution suggests a form for the 
observed luminosity function of N(K) **10*. In NGC 8060 the observed value of 
7 Is 7 ■ 0.89 a 0.03. In NGC 7083 and NGC 8083, the number of stars Is smaller 
and the uncertainties therefore larger, but consistent results are found of y - 
0.84 ± 0.05 and y * 0.81 ± 0.09 respectively. 

The relationship between the observed luminosity function. N(K), and the 
intrinsic luminosity function n(/f 0 ). where IC 0 is the absolute K magnitude, and 
ti(Kq) Is the number of stars In a magnitude interval dJTo, depends on a number 
of characteristics of the cluster, such as the spatial distribution of stars and 
dust. However, as long as n(/f 0 ) is everywhere the same in the nebula, then for 
N(K) ~ 10*. as was found for the three reflection nebulae. n(Ao) ~ 1 0 7jr °. This 
also assumes that any cutoff or change In slope in «(Xo) does not occur within 
the observed magnitude range of the sample from which N{K) was derived. For 
the finite number of stars observed here, it is possible that differences in the 
reddening of individual stars will alter the slope of the observed luminosity 
function. Since the reddening deduced from the color-color diagrams (Fig. 1 ), 
however, is typically only Ag ~ 0.5 mag, and at most Ag ** 1 mag, while the 
observed luminosity function is determined over 4 magnitudes at K, the 
reddening of individual stars should have little effect on the slope of the intrinsic 
luminosity function derived from the observations. 
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To convert the observed luminosity function n(Ao) Into a mass function one 
needs to know the relation between absolute K magnitude and stellar mass. 
This relation is completely unknown for the PMS stars, which comprise at least 
30% of the stars in NGC 2066. A simple-minded approach is to assume the mass- 
luminosity relation for PMS stars has the same slope as the mein sequence 
mass-luminosity relation. This assumption is more or less equivalent to 
assuming that the ratio of PlIS to main sequence stellar luminosity is 
approximately independent of spectral type within a given cluster. This 
simplistic assumption appears to be a fair approximation to published 2.2 fjm 
observations of PMS stars in other regions, for a wide range in spectral types, 
and to the H-R diagrams derived by Cohen and Kuhi (1079) for PMS stars, which 
show the PMS stars in a cluster to generally lie in a band on the H-R diagram 
parallel to and above the main sequence. 

A second problem in estimating the mass function is that the observed 
iuminosity function contains contributions from both PMS and main sequence 
stars. For this the simplifying assumption is made that the fraction of PMS stars 
is the same at each observed luminosity. While an increasing fraction of PMS 
stars with decreasing luminosity might be expected to be more appropriate, a 
constant fraction is consistent with the observations and more complicated 
expressions do not seem warranted. 

With these two simple-minded assumptions, then, the main sequence mass- 
luminosity relation combined with the derived cluster luminosity function can 
be used to estimate the cluster mass function. The main sequence mass* 
luminosity relation, for main sequence stars, is approximately Lg « M* *, where 
Lg is the luminosity at 2.2 fun and U is the mass. Since Lg « lO -04 * 0 , the 
luminosity function n(r 0 ) ~ 10 7jr ° implies a power law distribution with mass, 
n(Af) - M*. where fi is the power law index for mass, and n(Jtf) is the number of 
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stars in a logarithmic macs interval d logil. rrom tha obaarvad luminosity 
function in NGC 2060, with the above assumptions, one derives a mass power law 
index of fi * -1.8 ± 0.2. Ibis is not significantly different from the standard 
Salpeter mass power law index fi = -1.36, and agrees well with the range of mass 
power law indices observed in associations of T Tauri stars, fi = -1.35 to -2.9 
(Cohen and Kuhi 1979). 

d. 77ie attUar density and star formation affidtney 

Norman and Silk (1980) have described a scenario for star formation in 
dark clouds where stellar winds from T Tauri stars in the clouds provide the 
energy input for further formation of T Tauri stars. This theory explains the 
noncoeval star formation in dark clouds suggested by the observed spread of 
ages of T Tauri stars in the same cloud. Norman and Silk argue that densities of 
T Tauri stars above a critical density of 10 pc -8 are needed for this mechanism 
to work, and so it is of interest to calculate the density of T Tauri stars in NGC 
7023, 2023, and 2088. 

The spatial densities of the stars in NGC 7023, 2023, and 2088 can be 
estimated in two ways. First, one can assume spherical symmetry. Second, one 
can estimate the extent of the cloud along the line of sight by a comparison of 
the gas density n H and the gas column density N H . and then assume that the 
stars and gas are coextensive. The gas density and gas column density can be 
obtained from molecular line observations. 

For NGC 7023, the column density of 13 C0 is 1 x 10 1# cm* 1 (Loren. 
VandenBout, and Davis 1973), and the hydrogen density is ~ 10 s cm -8 (Whitcomb 
of of. 1981). For a 18 C0 to hydrogen ratio of 2 x 10 -8 (Dickman 1975), the 
approximate line of sight distance is ~ 2 pc. For NGC 2023, the column density 
of 18 C0 is 1 x 10 17 cm -8 (Milman sf of. 1975), and the hydrogen density implied 
by the C II observations (Pankonin and Walmsley 1978) is - 5 x 10* cm' 8 , so that 
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Um llna of tight distance it ~ 0.4 pc. For NGC 2066 the column density of “CO it 
6 x 10 1 * cm"* (Uilman tt al. 1075), and tha hydrogen dantlty la 4 x 10* cm"* 
(Loran 1961), implying a Una of tight diatanca of - 0.2 pc. Tha acan araaa in NGC 
7023, 2023, and 2066, for Um sample which hat been limited to at to be 
complete, are respectively 0.58, 1.4, and 1.3 pc*. and contain 6, 5, and 9 known 
PICS stars. Assuming spherical symmetry, the number densities of PUS stars 
are ~ 10, 3, and 6 pc"* in NGC 7023, 2023, and 2066. Adopting instead the line of 
tight distances from a comparison of N& and n* gives PUS stellar densities of ~ 
6, 10, and 30 pc"*, respectively. The actual densities of PUS stars are likely 
higher, since there are probably PUS stars fainter than the scan limits, and also 
the methods of this paper for identifying PUS stars provide only lower limits to 
their numbers. Thus, the densities of PUS start in these three regions are 
comparable to the densities required by the model of Norman and Silk (1960) to 
maintain continuing low mass star formations in dark clouds. 

A comparison of the column densities of stars and gas does not depend on 
uncertain estimates of the line of sight extent of the stellar cluster. An effective 
star formation efficiency is then the ratio of the mass per unit area in stars to 
the combined mass per unit area of stars plus gas. The most conservative 
estimate is to assume all stars whose spectral types are unknown are PUS stars, 
since these stars are less massive than a main sequence star of the same 
luminosity; to make no correction for possible reddening of the stars; and to 
include only stars already found in the scans limited so as to be complete. This 
estimate will therefore be a strong lower limit, especially since there are 
presumable fainter stars which contribute considerably to the mass of the 
cluster. The estimates of stellar mass contained in the areas scanned in NGC 
7023, 2023, and 2086 are respectively 10, 10, and 30 U Q under these 
assumptions. The masses are - 10-40% higher if stars in Table 2 below the scan 


completeness limits are Included. The gas mass per unit area is found from the 
“CO column densities above. The observed lower limits to the effective star 
formation efficiencies in NGC 7023, 2023, and 2368 are ~ 30%, 2%, ami 10% 
respectively. These efficiencies would be - 60%, 3%. and 20% if stars not known 
to be PlfS stars were assumed to lie cm the main sequence. While the conversion 
of gas into stars in NGC 7023 appears to be relatively efficient, this reflects in 
part the fact that the NGC 7023 stellar cluster is not at a molecular peak. In all 
three nebulae the molecular gas extends over a larger area than does the stellar 
cluster, so that the total star formation efficiency for the cloud is much lower 
than the local efficiency within the cluster. It would be extremely interesting to 
examine how the efficiency of star formation varied with radius in the stellar 
cluster. This comparison, however, must await molecular line observations with 
higher resolution than the angular resolution of ~ 3' presently available. 


Tha properties of atari found in 2.2 pm warn of redaction nabuiaa hava 

boon investigated. and several oonclualona about the nature of the atara can be 

drawn: 

1 . The number of atara found la much larger than the number of expected 
background -tars, Implying that almoat all of the atara found are in fact 
aaaoclated with the cloud. 

2. Consideration of a number of indicators of the pre-main sequence (PMS) 
nature of stars, such as hydrogen line emission, infrared excesses, or 
Irregular variability, shows that at least 30-60 X of the stars found in these 
stellar clusters have not reached the main sequence. 

3. The observed variability Indicates that at least 50 X of the stars in NGC 7023 
are PMS objects, but most of the nebular variable stars in NGC 7023 do not 
have Identifiable infrared excesses. This suggests that identification of 
young stars by their infrared excesses, as was done in NGC 2023 and 2068, 
severely underestimates the true numbers of PMS stars. 

4. Tha observed luminosity function of the young stellar clusters is found to be 
N{K) - 10**. with 7 « 0.29 £ 0.03. 

5. The star formation efficiency local to the stellar cluster is found to be at 
least 30 X in NGC 7023, although the overall efficiency for the entire dark 
cloud is lower. 

6 . The spatial densities of PMS stars are found to be 2 3*30 pc' 8 , comparable 
to values required by the model of Norman and Silk (i960) to maintain 
constant low mass star formation in the dark clouds. 
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7. The surface density of the stars In all three nebulae decreases as the 
inverse of the projected radius from the cluster center. This, if spherically 
symmetric, corresponds to a volume density which depends on the true 
radius as r~*. 
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Table 1 


$• (pacing 4 «ean 

n n difction* 


A scan limits 

ium) a* g* mMm* data 

a. Scans of NGC 7023 


1.65 

23E-23W 

23E-23W 

23S-158S 

23N-203N 

11.3 1979AUG22 

11 

20 

r& 

2.2 

158E-113W 

113W-158W 

203N-203S 

23N-113S 

11.0 1979AUG23-24 

11 

20 

ra 

2.2 

113W-158W 

113W-156W 

23E-6BE 

23N-203N 

113S-203S 

23S-68N 

11.0 1979SEP03 

11 

32 

ra 

2.2 

158E-203E 

203W-158W 

203E-316E 

203W-312W 

203S-203N 

203S-U3N 

68S-23N 

68S-23N 

11.7 1981JUL30»AUG01 

11 

23 

ra 

2.2 

129E-424E 

462W-301W 

81W-81E 

14BE-1 7 2E 

15N-176N 

157S-4N 

01S-81N 

214S-52S 

12.0 1981N0V15-21 

16 

129 

dac 




6. Scans of NGC 2023 




2.2 

250E-250W 

273N-200S 

11.1 19B1JAN05-08 

20 

70 

dec 

2.2 

50W-50E 

100S-100N 

11.1 1981JAN07 

10 

70 

dec 

2.2 

136E-310E 

228W-307W 

ie3S-6BN 

23S-68N 

11.0 1981FEB01 

11 

36 

ra 

2.2 

/ 

/ 

/ 19810CT20 

19 

344 

ra 




e. Scan* of NGC 2066 




1.65 

0W-382W 

382S-67N 

11.0 1961FEB03-05 

19 

344 

ra 

2.2 

0W-382W 

3B2W-764W 

29N-191N 

48S-191N 

10.6 1981FEB06-12 

19 

344 

ra 

2.2 

50E-50W 

0N-100N 

12.6 1961SEP16-17 

10 

60 

dac 

2.2 

50E-250E 

ON-lOON 

11.9 1981SEP17-19 

20 

60 

dac 

2.2 

0W-3B2W 

0S-1B1S 

8.7 19810CT13-14 

19 

344 

ra 

2.2 

90E-0E 

100S-200S 

11.9 1981N0VU 

20 

55 

dac 

2.2 

161E-161W 

0W-322W 

322E-161E 

181E-81E 

OS-322S 

161N-0N 

0S-161S 

205N-44N 

12.0 1981N0V17-20 

16 

129 

dac 


(a) Offset in arc see from HD 800779 * star A In NGC 7083, from HD 37003 * star 
A in NGC 8083, and from HD 38863-N ■ star 33 In NGC 8066. 

(b) Faintest magnitude detectable 

(c) Diameter of diaphragm 

(d) Spacing between source and sky positions 

(e) The scan technique used was a raster scan covering the area indicated. The 
direction of the scan lines was along the line connecting the source and sky 
positions. 

(f) Specific search for a particular star rather than a systematic scan. 





found m 2.2 pm scant of NGC 2023 



found vn 2 2 /an and 1 . 06 fim Meant of NQC 2066 
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TO*«3a 


Variability in NGC 7023 


Name 

K* a( XT 


date 0 — 

36 

9.50 


1972 

36 

11.0 


1981JUL31 

36 

11.0 


1981AUG01 

36 

11.0 


1981AUG01 

36 

11.05 

11.51 

1981NOV15 

J 

9.2 


1979AUG23 

J 

9.4 


1979AUG24 

J 

9.3 


1979AUG24 

J 

9.58 

10.08 

1979SEP02 

J 

8.97 

9.62 

1981JUL27 

J 

9.1 


1981NOV20 

K 


10.0 

1979AUG22 

K 

9.4 


1979AUG23 

K 

10.2 


1979AUG24 

K 

9.40 

9.84 

1979SEP02 

K 

9.36 0.06 

9.80 

1981JUL27 

K 

9.6 


1981NOV19 

E 

10.50 


1972 

E 


11.3 

1979AUG22 

E 

10.5 


1979AUG23 

E 

10.7 


1979AUG24 

E 

10.60 

11.13 

1979SEP02 

E 

10.68 

11.30 

1981JUL26 

E 

11.4 


1981NOV19 

L 


10.1 

1979AUG22 

L 

9.4 


1979AUG23 

L 

9.5 


1979AUG23 

L 

9.6 


1979AUG23 

L 

9.6 


1979AUG24 

L 

9.55 

10.05 

1979SEP02 

L 

9.30 

9.64 

1981JUL27 


-32- 


Tabie 3* cdotiauad 

Variability in NGC 7023 


a m> ■ s. tuoi «• mi d ‘ i ‘‘ 


M 

11.0 


1979AUG23 

M 

11.0 


1979AUG24 

M 

10.99 

11.29 

1979SEP02 

M 


11.42 

1981JUL29 

M 

11.01 

11.37 

1981JUL30 

M 

11.8 


1981N0V19 

D 

10.2 


1979AUG23 

D 

10.04 

10.30 

1979SEP02 

D 

10.2 


1979SEP03 

D 

10.38 

10.58 

1981JUL28 

D 

10.53 

10.68 

1881JUL26 

D 

10.3 


1981N0V19 

C 

10.50 


1972 

C 

9.9 


1979AUG23 

C 

10.1 


1979AUG24 

C 

10.1 


1979AUG24 

C 

10. 0B 

10.88 

1979SEP02 

C 

10.25 

10.84 

1981JUL28 

C 

10.4 


19B1N0V19 

T 

11.19 

11.49 

1979SEP03 

T 

11.23 

11.47 

1981JUL29 

T 

11.0 


1981JUL30 

T 

11.9 


1981N0V18 


(a) K and H magnitudes given to 0.1 mag are scan magnitudes, and are 
uncertain by £ 0.3 mag. K and H magnitudes given to 0.01 mag are 
photometric magnitudes. 

(b) Uncertainties are ± la, and include both photometric and statistical 
uncertainties. The uncertainties are given only when greater than 0.05 mag 
for photometry; no uncertainties are given for scan magnitudes. 

(c) Dates for which the year only is given are for observations by other 
observers, and refer to the year in which the data were published. The 
supplemental observations were taken from Strom et al. (1972). 
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Trtteab 


KsrLabiLity in NGC2023 


F 

8.65 0.28 

1975 

F 

8.5 

19B1JAN06 

F 

8.7 

1981 JAN 06 

F 

9.9 

1981FEB01 

F 

9.28 

1901SEP15 


(a) K magnitudes given to 0.1 mag are scan magnitudes, and are uncertain by 
2 0.3 mag. K magnitudes given to 0.01 mag are photometric magnitudes. 

(b) Uncertainties are t lo, and include both photometric and statistical 
uncertainties. The uncertainties are given only when greater than 0.05 mag 
for photometry; no uncertainties are given for scan magnitudes. 

(c) Dates for which the year only is given are for observations by other 
observers, and refer to the year in which the data were published. The 
supplemental observations were taken from Strom et a l. (1975), and are 
the sum of the contributions of S102, S103. S104, and 5110 (see footnote to 
Table 26). 
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1febte3e 


Variability in NGC2068 


Name A« o(K\' H* a(H)' 


4 

10.50 

0.25 



4 

9.09 

0.1 

10.25 

0.1 

4 



10.6 


11 

9.B 




11 

11.1 




13 

10.60 

0.25 

11.54. 

0.25 

13 

10.2 




13 

11.0 




14 

9.40 

0.2 

9.70 

0.2 

14 

9.91 

0.1 

10.59 

0.25 

14 

10.0 




15 

9.3 




15 

11.91 

0.09 

12.75 

0.09 

15 

12.0 




16 

9.7 




16 

11.07 

0.25 



16 

11.2 




16 

11.2 




17 

9.0 




17 

10.08 

0.1 

11.04 

0.25 

17 

9.9 




17 

9.6 




17 

9.9 




17 

9.8 




17 

9.63 


10.40 


17 

9.8 




17 

9.78 


10.49 


19 

9.2 


10.3 


19 

9.94 

0.1 

10.35 

0.1 

19 

9.4 




19 

9.9 




19 

9.0 




19 

9.3 




19 

10.1 




19 

9.80 


11.23 


19 

9.67 


11.06 



1975 

1978 

1981FEB04 

1975 

1981NOV20 

1976 

19B1FEB12 

1981NOV20 

1975 

1978 

1981NOV20 

1981FEB12 

19B1N0V19 

19B1NOV20 

1975 

1978 

19B1N0V1B 

1981N0V18 

1975 

1978 

1981FEBG6 

19B1FEB12 

19B1FEB12 

19B1N0V18 

1981NOV19 

19B1NOV20 

1982JAN18 

1975 

1976 

1981FEB06 
19B1FEB06 
19B1FEB12 
1961 FEB 12 
1981N0V1B 
1981N0V19 
1982JAN16 
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Table 3c continued 
WtriaJbility in NGC8068 


a. rtm* 


20 

8.90 

0.2 

20 

9.64 

0.1 

20 



20 

9.9 


20 

9.7 


26 

B.B 


26 

9.87 

0.1 

26 



26 



26 

9.6 


26 

9.6 


26 

10.1 


26 

9.74 


26 

9.B6 


26a 

9.60 

0.2 

26a 

10.1 


26a 

10.15 


26b 

10.00 

0.2 

26b 

10.5 


26b 

10.62 


40 

9.5 


40 

10.1 


40 

10.1 


40 

10.7 


47 

9.5 


47 

10.33 

0.1 

47 



47 

10.4 


47 

9.9 


47 

10.6 


47 

10.2 


47 

10.3 


47 

10.4 



MLmm 

qua* 


10.20 

0.1 

1975 

1976 

10.6 


1981FEB03 

10.10 

0.1 

19B1N0V18 

19B1NOV20 

1975 

1976 

10.6 


1981FEB03 

10.3 


1981FEB03 

10.15 


1961FEB12 

1981FEB12 

1981N0V17 

19B1NOV19 

10.48 


19B1NOV19 

1975 

1981NOV18 

1982JAN16 

11.00 


1975 

1981NOV18 

1982JAN18 

10.B8 

0.25 

1975 

1961SEP16 

1981NOV11 

1981NOV17 

1975 

1976 

11.0 


1981FEB03 



19B1FEB12 

1981FEB12 

19B1SEP17 

1981SEP19 

1981NOV17 

1981NOV20 


(a) K and H magnitudes given to 0.1 mag are scan magnitudes, and are 
uncertain by 2 0.3 mag. K and H magnitudes given to 0.01 mag are 
photometric magnitudes. 

(b) Uncertainties are ± la, and include both photometric and statistical 
uncertainties. The uncertainties are given only when greater than 0.05 mag 
(or photometry; no uncertainties are given (or scan magnitudes. 

(c) Dates (or which the year only is given are for observations by other 
observers, and refer to the year in which the data were published. The 
supplemental observations were taken from Strom at at. (1975), and 
Strom, Strom, and Vrba (1976). 


Figure Captions 


Fig. 1. J-H v*. H—K colors in magnitudes of ths stars observed in (a) NGC 7023, 
(6) NGC 2023, and (c) NGC 206B. Error bars are ± Iff. Tilled circles indicate 
observations from this paper, while open circles in Tig. 2(c) indicate 
observations from Strom. Strom, and Vrba (1976). The solid curve is the m ai n 
sequence, from AOY to M5V. The solid line indicates the expected reddening line 
for a background K giant, while the dashed line indicates the reddening line for 
an AOY star. 

Tig. 2. Plot of positions of stars found in (a) NGC 7023, (b) NGC 2023, and (c) 
NGC 2068. Pluses indicate stars found at 2.2 or 1.65 /xm, open circles indicate 
stars found at 0.655 or 0.8 fi m, an asterisk indicates the visual illuminating 

star of each reflection nebula. The solid outline indicates the scan area for the 
scans at 2.2 /xm from which complete samples of stars were drawn, while the 
dashed outline indicates additional area scanned at 2.2 /xm. 

Tig. 3. Cumulative counts N{K) vs. observed K magnitude for the stars in NGC 
2068. The cumulative counts are given as number of stars per square degree 
brighter than the observed K magnitude. Error bars are ± Iff. The indicated 
line is a least squares fit of the function N(K) = No 10 Tjr to the data. 
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Figure 2a 
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Figure 2c 
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